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Abstract 

o, 

The eigen mode of spin oscillations with ui ~ ^58/35 A is predicted to exist besides 
already known spin waves with oj ~ A/y5 in the triplet superfluid neutron con- 
densate in the inner core of neutron stars. The new mode is kinematically able to 
decay into neutrino pairs through neutral weak currents. The problem is considered 
in BCS approximation for the case of 3 i-2 — 3 -^2 pairing with a projection of the total 
angular momentum m,j = which is conventionally considered as preferable one at 



supernuclear densities. 
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A superfluidity of the inner core of neutron stars plays a crucial role in theirs 
cooling scenario. The energy gap A arising in the quasiparticle spectrum be- 
low the critical condensation temperature T c suppresses the most of neutrino 
emission mechanisms [lj . According to the minimal cooling paradigm [2|3Pf5] , 
under these conditions, the most efficient energy losses from the star volume 
can take place at a recombination of thermal excitations in the form of broken 
Cooper pairs. The neutrino emission at the pair-recombination processes oc- 
curs through neutral weak currents in the axial channel of weak interaction^ 
and can be very efficient, in the triplet superfluid neutron liquid, a few below 
the critical temperature |S]. However, the corresponding neutrino emissivity 
falls rapidly with lowering of the temperature because the number of broken 
pairs, having the excitation energy larger than 2 A, decreases exponentially. 
In this case the collective excitations of the condensate can dominate in the 
neutrino production. 

Since we assume that the condensate consist of neutron pairs in the triplet 



1 The vector channel of the neutrino radiation through neutral weak currents is 
strongly suppressed in the non-relativistic case [7|8j . 



Preprint submitted to Elsevier 21 July 2011 



state it is natural to expect the collective modes associated with spin fluctu- 
ations of the condensate 2 1. Such collective excitations with the energy lower 
than 2A might undergo the weak decay into neutrino pairs. Recently spin 
waves with the excitation energy u = A/v5 was predicted to exist in the su- 
perfluid spin-triplet condensate of neutrons [8f9][T0] . Because of a rather small 
excitation energy, the weak decay of such waves leads to a substantial neutrino 
emission at the lowest temperatures T <C T c , when all other mechanisms of 
the neutrino energy losses are killed by the superfluidity. 

In Refs. [8.9.10J, the eigen-mode of spin oscillations in the 3 P 2 superfluid neu- 
tron liquid was studied in a simple model restricted to excitations of the 
condensate with I = 1. In this paper we demonstrate that extending of the 
decomposition up to I = 1, 3 leads to a very small frequency shift of the known 
mode, u = A/\/5, but opens the new additional mode of spin oscillations with 
the finite energy gap u> (q = 0) < 2A. The problem is considered for the case 
of 3 P2— 3 F 2 pairing with a projection of the total angular momentum rrij = 
which is conventionally considered as preferable one at supernuclear densities. 

We will examine the spin modes within the BCS approximation 3 1. Let us 
remind briefly the theory of spin density excitations in the condensate. The 
order parameter, D = D a/3 , arising due to triplet pairing of quasiparticles, 
represents a 2 x 2 symmetric matrix in spin space, (a, =j, i)- The spin-orbit 
interaction among quasiparticles is known to dominate in the nucleon matter 
of a high density. Therefore it is conventional to represent the triplet order 
parameter of the system as a superposition of standard spin-angle functions 
of the total angular momentum (j,rrij), 

&,l,mi) f^ _ V- (\\. 
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Assuming that the pair condensation occurs into the state with a total angular 
momentum j = 2 we use the vector notation which involves a set of mutually 
orthogonal complex vectors bj m . (n) defined as 



■j 



k mj (n) = ~Tr {g*&> 1 ™) , b,,_™, = (-P K mj , (2) 

where & = (<Ti,<72> 03) are Pauli spin matrices, g = 1&2, and the angular 
dependence of the order parameter is represented by the unit vector n = p/p 



2 Previously spin modes have been thoroughly studied in the p-wave superfluid liq- 
uid 3 He with a central interaction between quasiparticles |ll|12|13|14|15j . These 
results cannot be applied directly to the triplet-spin neutron superfluid condensate, 
where the pairing is caused mostly by the spin-orbit interaction between quasipar- 
ticles (see details in Ref. [9]) 

3 Throughout this paper, we use the system of units H = c = 1 and the Boltzmann 
constant ks = 1. 



which defines the polar angles (0,<p) on the Fermi surface. The vectors b/ >m . 
are mutually orthogonal and are normalized by the condition 



\b*',r^b/,m.;) - Su'S m . m '.. (3) 

Hereafter the angle brackets denote angle averages, (...) = (47r)~ J dn.... 

The block of interaction diagrams irreducible in the channel of two quasipar- 
ticles, r Q/3/y< 5, is usually generated by expansion over spin-angle functions. The 
spin-orbit interaction among quasiparticles is known to dominate at high den- 
sities. This implies that the spin s and orbital momentum 1 of the pair cease 
to be conserved separately, and the complete list of channels includes the pair 
states with j = 0, 1,2, and \rrij\ < j. These nine complex states exhaust the 
number of independent components in the matrix order parameter arising at 
the P-wave pairing caused by the strong spin-orbit forces. The pairing in the 
j = 2 channel dominates, and due to relatively small tensor components of the 
neutron-neutron interaction the condensation of pairs occurs in the 3 P2+ 3 P2 
state. In this pairing model, contributions from 3 P2 - > 3 Pq or 3 P2 — > 3 P\ transi- 
tions are deemed to be unimportant. Such assumption is somewhat vulnerable 
especially when considering excited state of the condensate. Unfortunately the 
detailed information on the in-medium effective interaction between neutrons 
in the channels j = 0, 1 is currently unavailable and requires a special inves- 
tigation. Hence we take the approximation to neglect the j = 0, 1 coupling 
throughout this paper. From now on we omit the suffix j everywhere by as- 
suming that the interaction occurs in the state with j = 2. Thus we assume 
I — j ± 1, and 

gT a ^ 5 (p, p') = £ (-l)^V w (p,p')(b lmj (n)ag)^(g&b* llm .(n'))^, (4) 

I'lmj 

where Vw ij>,p') are the interaction amplitudes, and 1,1' = 1,3, in the case of 
tensor forces; g = pfM*/tt 2 is the density of states near the Fermi surface 
in the normal state. The effective mass of a neutron quasiparticle is defined 
as M* = pf/vf, where vf ^ 1 is the Fermi velocity of the non-relativistic 
neutrons. 

The order parameter is of the following general form 

D(n)=J2^i,m J (*b Wmj )g. (5) 



The ground state occurring in neutron matter has a relatively simple structure 
(unitary triplet) [TBfTT] . where 

^A l)m; b []m] (n) = Ab(n). (6) 



On the Fermi surface, A is a complex constant, and b (n) is a real vector 
which we normalize by the condition 



;& 2 (n)) = l. (7) 

The following orthogonality relations are also valid: 



J WZ^tmj^l',^'- - Smjmfilmhl's 



J 7^7 ( bb *,"J ( bb i',m$) = 6 rn 3 mr ( b K, mj ) ( bb F, mj ) • (9) 

Thus the triplet order parameter can be written as 

D (n) = Abfrg. (10) 



Making use of the adopted graphical notation for the ordinary and anomalous 
propagators, G = — » — , G~(p) = — « — , F^ = » « , and F^ = « » , it 
is convenient to employ the Matsubara calculation technique for the system in 
thermal equilibrium. Then the analytic form of the propagators is as follows 
[T8HT91 



G {pr,, p) = G (p v , p) <W , G (p rj ,p) = G (p v ,p)5 a p, 

F™ (p„, p) = F (p v , p) bag , F® (^ p ) = F (p v , p) gab, (11) 

where the scalar Green's functions are of the form G~ (p v , p) = G (— p v , — p) 
and 

ofa - rt = ^^- F ^ = g^- (12) 

Here p v = i7r(2r)+l)T with rj = 0, ±1,±2... is the Matsubara's fermion 
frequency, and e p = p 2 / (2M*) —p 2 F / (2M*). The quasiparticle energy is given 
by E 2 = e 2 +A 2 b 2 (n), where the (temperature-dependent) energy gap, Ab (n), 
is anisotropic. In the absence of external fields, the gap amplitude A is real. 

Finally we introduce the following notation used below. We designate as 
Ixx' {oj, n, q; T) the analytical continuations onto the upper-half plane of com- 
plex variable u of the following Matsubara sums: 

l X x< (w«, n, q; T) = T J^ - J de p X (p v + u K ,p+-J X' (p v ,p--J , (13) 

where X, X' e G,F,G~, and u K = 2mTn with k = 0, ±1, ±2. ...These are 
functions of u, q and the direction of the quasiparticle momentum p = pn. 

We will focus on the processes with u> 2 < 2A 2 b 2 and with a time-like mo- 
mentum transfer, q 2 < u 2 . In this case the key role in the response the- 
ory belongs to the loop integral Iff- A straightforward calculation yields 



Iff ( n , w, Q n ; T") = X (n,w; T) + O {q 2 v 2 F /u 2 ), where 

f°° de A 2 E 

MW,T)=l e UP - {u + «? ***&■ <14) 

Insofar as q 2 v 2 F /u 2 <C 1 and q 2 v 2 F /A 2 C 1 we will neglect everywhere small 
corrections caused by a finite value of space momentum q. 




*,** (p) = - £ ^ / tilT^Vu' (p, J/) a bo /t E 



Pl + El, /• 

(15) 
We are interested in the processes occurring in a vicinity of the Fermi surface. 
To get rid of the integration over the regions far from the Fermi surface we 
renormalize the interaction as suggested in Refs. [251126] : we define 

VP (p,p'; T) = V lv (p,p f ) - W ^f-V ip,p") (GG-)" VS (pV; T) , 

I" 

(16) 
where the loop (GG ) n is evaluated in the normal (non-superfluid) state. In 
terms of V^, the renormalized gap equations can be written in the following 
matrix form 

Am,A_ A (vi; > -^'U(bb* m ^;, 



VJ V-^13 ^33 / \\ uu 3,m 
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assuming that in the narrow vicinity of the Fermi surface the smooth functions 
Viii (p,p ! ) and A (p') may be replaced with constants. In obtaining Eq. ( TTTj) 
the fact is used that the interaction matrix is symmetric on the Fermi surface, 
V31 = V13. The function A (n) arises due to the renormalization procedure. It 
is given by 



A (n) = I r dc I \ tailh ^±5? _ I tanh 1: 

K ' 2 Jo We 2 + A 2 b 2 2T e 2T . 



The interaction matrix can be diagonalized by unitary transformations V 
t/W^with U being an unitary matrix 



U = (U-'y = -—-r \ V _ V _|, (19) 




where V ± = J (v&> - V^ + AV^ 2 ± (vg - V £>). 



One has UVW = diag (W-, W+) with 



(20) 



Applying the unitary transformation U to the gap equations ( TT7I) yields two 
coupled equations: 



V+A hm . + JV-A 



-w-£a 



Lm,i 



^3,771, 



V + b* t b,, ro . + x/VLb 



3,m,j "!,m- I ^ 



(21) 



-^+E A 



V^A 



3,m, 



(,m,- 



^-bl/m.b^m, - VV+bg b. 



'3,mj °l,m,j 



A 



(22) 



In obtaining these equations we made use of Eq. fl6]) and orthogonality relations 
OH]), assuming that the energy gap b (n) A is azimuth-symmetric |16fl7f20f21|l22|l23p4 



We are interested in the linear medium response onto the external axial- vector 
field. The field interaction with a superfluid should be described with the 
aid of two ordinary and two anomalous three-point effective vertices. In the 
BCS approximation, the ordinary axial- vector vertices of a particle and a hole 
are to be taken as a and <x , respectively. The anomalous effective vertices, 
T^ (n;w, q) and T^ 2 ^ (n;w, q) are given by the infinite sums of the diagrams 
taking account of the pairing interaction in the ladder approximation |27j . 
These 2x2 vector matrices are to satisfy the Dyson's equations symbolically 
depicted by graphs in Fig. [TJ Analytic form of the above diagrams is derived in 
Refs. p]. After some algebraic manipulations the BCS equations for anomalous 
vertices can be found in the following form (for brevity we omit the dependence 
of functions on u and q) : 



f « (n) = 5>Wn)$5>«'2 ( X ^- Tr [9 (*b?, m .) T (1) 



Irrii 



— XppTr 



<xb* m J («xb) gf < 2 > (<xb)] - ^-l FF 2i (b?, m . xb) V (23) 



f ( 2 ) (n) = E^K m .(n) E V„,- (x G - G Tr [(<xb, m .) gf 

Irrij I' 



<2) 



-1 ff Tt 



&b l/m .)(ab)T^g(ab 



UJ 

A 



I FF 2i [bii m 



l m . Xb 



(24) 
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Fig. 1. Dyson's equations for the anomalous vertices. The shaded rectangle repre- 
sents the pairing interaction. 

Inspection of the equations reveals that the anomalous axial-vector vertices 
can be found in the following form 

f « (n, w) = J2 B ^ h (w) (*b, lt J & (25) 



T( 2 )(n, W )=£B ; , m >)s(*b* mj ). (26) 

As explained above we are interested in solutions with q = 0. Then inserting 
of these forms into Eqs. (1231) . (1241) allows to obtain the equations for B/ jmj (u>). 
We write the result in the matrix form (For brevity we omit the dependence 
on n and cu) 



Bi, 
Bs 



'j 






Ei(2b(Bb^)(Sb i|ny ))B 



-2| * : "\ ;-\' v '*':' |+|i 



/ 



^obL.xb 




E* (X (bb^.) (bb,, m .) ) B iim . j A \(lo (b^ m . xb 



In this equation, the interaction matrix can be diagonalized by the unitary 
transformation f|T9|) . Further simplification is possible due to the fact that 
by virtue of Eqs. ( 12T1) . (|2"2"|) the coupling constants PV± can be removed out 
of the equations. Explicit evaluation of equations obtained in this way for 
arbitrary values of u and T requires numerical computation. However, we 
can get a clear idea of the behavior of the vertex functions using the angle- 
averaged energy gap A 2 b 2 — > (A 2 b 2 ) = A 2 in the quasiparticle energy E p . In 
this approximation, the functions X (u, T) and A (T) can be moved beyond 
the angle integrals. Performing trivial integrations we then get a set of linear 
equations (two equations for each value of rrij). It is convenient to denote 

ffi = ((bi im ,6) (K^b)) , (28) 



and 



2A 



Then the set of equations can be written in the form 



Bi >mj . 









-in V^ <bU x b + V F - b 3,-, x b 



(29) 



(30) 



Bi >m 
+B3™ 



(™i) 



-^(n 2 -/^ 



j( m j) 



v + iw-pw') + jv-p% 



v+p^ 



jK) 



-in - jy_ b* x b + jy + b* . x b 



which can be solved to give 



Bi, 



-in r 



B3, m . 



X 

-in r 



(fi 2 -^)(b;, ra] xb) + tf(b; im] xb 



X 



^2 «W( /^ 



n^ - py? 



3,m, 



X b) + ^ } (Km* X b 



with 



x (n) ^ n 4 - n 2 (p™ + &) + ^ ] & - ffipft 



(31) 

(32) 
(33) 
(34) 



As is well known, poles of the vertex function correspond to collective eigen- 
modes of the system. Eigen- frequencies, n = VlS m3 \ of such oscillations satisfy 
the equation x (n^ m ^ j = 0. This equation gives 

(ni mj) ) 2 = \ (pti ] + P& ± \l{Pt i) -&) 2 + *&P&^ ■ (35) 

Notice that the interaction parameters, V±, drop out of the above solutions, 
which depend explicitly only on the partial gap amplitudes. This means that 
the contribution of excited bound pairs with / = 3 into the spin oscillations is 
caused basically by spin-orbit interactions but not by the tensor forces. 

Indeed, in Eqs. ( 152]) - (j55|) . the equilibrium order parameter is specified solely 
by means of the real vector b. If we switch off the interaction in the 3 F2 and 
3 P2— 3 F2 channels and consider pure 3 P2 pairing with rrij = we are then left 
with b = bio and A = Ai . In this case, in Eqs. ( 132]) . (133]) . one has: 



/£fc 



xb 



/ % 0U x b) = for m, = 0, ±2, (36) 



and the non-trivial solutions exist only for rrij = ±1. The explicit form of b^ m . 
can be obtained from Eq. fl2]): 



bi, 



/ \ 

-m 

-n 2 
2n 3 



I 



bi. 



y -^3 j 



n 3 

in 3 
ni + in>2 



(37) 



b 3 ,o 



m (1 - 5n 2 3 ] ^ 



n 2 (1 - 5n|) 
\n 3 (3-5nl) ) 

Making use of these expressions in Eq. 



, b 3 i 



n 3 (1 - 5ni (ni + in 2 )) 
in 3 (1 + 5in 2 {n x + in 2 )) 

{m + in 2 ) (1 
we find 



5n|) 



(3J 



Pl,l " go ' ^3,3 - ^ 



Pl,3 



P3,l 




(39) 



Inserting these values into Eq. $2} we find 4/?f 3 1) /3f 1 1) < (/3f 3 1} - Mf^) 2 . By 

neglecting the small term Af3 1 3 (3 3 x under the root in Eq. fl35|) we obtain two 
(twofold) eigen-frequencies of spin oscillations in the condensate with rrij = 
±1: 



ijj 



(rrij) 



2Adf3[f ) = 4=A 



UJ 



(rrij) 



2AJ/3 



? (±i) 

; 3,3 




(40) 
(41) 



In Refs. [9][T0] . eigen- modes of spin oscillations in the 3 P 2 superfluid neutron 
liquid was studied in a simple model restricted to excitations of the condensate 
with 1 = 1. The spin wave energy (at q = 0) was found to be u> mj = A/y^- 
Equations ( |35|) . ( 140]) . (|4T|) show that extending of the decomposition up to 
/ = 1, 3 in Eqs. (1251) . (126|) leads to a very small frequency shift of the known 



mode, u 
with uj = uj 



(rrij) 



id 
(rrij) 



UJ 



rrij, 



but opens the new additional mode of spin oscillations 



Neutrino decays of spin waves can play an important role in the cooling sce- 
nario of neutron stars. A simple estimate made in Ref. [10] has shown that 
the decays of spin waves with co m . = A/v^5 can become the dominant cooling 
mechanism in a wide range of low temperatures and modify the cooling tra- 
jectory of neutron stars. As well as the first mode, the second mode of spin 
oscillations is kinematically able to decay into neutrino pair. Therefore let us 
examine the wave excitation energies more accurately with taking into account 
the tensor forces. We will again focus on the condensation with rrij = by 



assuming A 2 = A 2 + A§ , and 



b(n) 



Ai, 
A 



b i,o (n) 



A 



3,0 



A 



b 3 ,o (n) 



(42) 



In this case Eqs. (1361) are still valid and the non-trivial solutions to Eqs. 
( 1331) exist only for rrij = ±1. Insertion of the expression (1421) into Eq. 
results in 



(±1) _ 

Pl.l - nn a 2 



j(±l) 
J 3,3 



P3,3 — ™ A 9 




& 



/#e 



(1) 



7^AT 



58Ajf 

7A? 



32 ^A 3 28 Af 



1 /R A * 



^A 3 32 A 2 / 

i-^aT + Ta!, 



Because /3^, = /? 



U' 



/3; F we further omit the superscript rrij 



assuming that all the frequencies are twofold. Making use of Eqs. 

' 1S iV6 As 



u 



K 



14 



23 A 2 

'Ax + 21 A 2 



15 25 ^A 3 485 Al 
^28~49 V6 aT + 147 A? 



370 /-Al 55 Af 
441 V Af + 63Af 



(43) 

±1 by 
we find 



(44) 



u. 



A 2 



13 
14 



1 ^A 3 23 A^ 



^A 2 

s^AT + ilA! 



+, 



15 25 ^A 3 485 A^ 



-— v 7 ^ 



+ 



370 ^ A 



+ 



55 At 



\ 28 49 Ai 147 A 2 441 A? 63 Af 



(45) 



In Fig. 121 the energy of the collective spin excitations (at q = 0) is shown vs 
the ratio of the partial gaps x = A 3j o/Ai ;0 . According calculations of different 
authors, at the Fermi surface one has A 3 ~ 0.17Ai (see, e.g., Ref. |20j). In this 
case our theoretical analysis predicts two degenerate modes with to = uj- = 
0.42 A, and two degenerate modes with u = u + = 1. 19 A. 

Because of a rather small excitation energy the decay of the corresponding col- 
lective spin excitations into neutrino pairs should lead to an extension of the 
low-temperature domain where the volume neutrino emission dominates the 
surface gamma radiation in the star cooling. This effect was already demon- 
strated in Ref. [10] , where only the lowest branch of the collective spin excita- 
tions u = A/v/5 — uj- has been taken into account. The neutrino emissivity 
caused by the decay of the new spin modes predicted in this paper will be 
considered elsewhere. 
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0.3 0.4 

'3,0^1,0 



X=Aon/A- 



Fig. 2. The energy gaps for the collective spin excitations u_ J and uj + J vs the 
ratio of partial gap amplitudes in the 3 i*2 and 3 i-2 channels. The energy gap of a 
neutron quasiparticle is given by A 2 = Af + A| . 
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